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Abstract 

A novel view for the emergence of chaos in Lorenz-hke systems is presented. 
For such purpose, the Lorenz problem is reformulated in a classical mechan- 
ical form and it turns out to be equivalent to the problem of a damped and 
forced one dimensional motion of a particle in a two-well potential, with a 
forcing term depending on the "memory" of the particle past motion. The 
dynamics of the original Lorenz system in the new particle phase space can 
then be rewritten in terms of an one-dimensional first-exit-time problem. The 
emergence of chaos turns out to be due to the discontinuous solutions of the 
transcendental equation ruling the time for the particle to cross the inter- 
mediate potential wall. The whole problem is tackled analytically deriving a 
piecewise linearized Lorenz-like system which preserves all the essential prop- 
erties of the original model. 

PACS number: 05.45.-a 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The Lorenz dynamical system, originally introduced by Lorenz in order to describe 
in a very simplified way the Rayleigh-Benard problem [0,0, immediately became important 
in itself as one of the most studied low-dimensional chaotic systems. Still today the Lorenz 
model represents a paradigmatic example for both theoretical and numerical investigations 
in checking some results in chaos theory 0-0, in the study of geometrical properties of 
dynamical systems P- p!l|] , in nonlinear analysis of time series [|1^,0], in the stabilization 
and synchronization of coupled systems [|l^-[l6[|, and so on. 



Nevertheless, despite the great attention attracted over past decades, some fundamental and 
rigorous results have been obtained quite recently, as for instance the proof of the existence 



of the Lorenz attractor [^,|T^], usually using somewhat sophisticated mathematical tools. 

On the contrary, our aim here is to provide a description of Lorenz system dynamical 
features, which requires quite simple analytical tools and, at the same time, allows a very 
intuitive inspection in Lorenz-like chaos. Preliminary results have been reported in a short 



communication |jT9l. Here we shall give more details and new results. 



Our interpretation will base itself upon the fact that in the "steady state", i.e., far from 
the initial transient and when the memory of the initial conditions has been lost, the Lorenz 
system is equivalent to a suitably constructed second order integral-differential equation. 
This equation can be regarded, for instance, as a customary second order one-dimensional 
classical mechanics equation with a peculiar forcing term. The corresponding dynamics can 
be interpreted as the one- dimensional motion of a particle in a conservative quartic two- well 
potential, subjected to a viscous damping and to an additional force resulting from the past 
history of the motion. The latter force turns out to be essential for chaos to emerge as it 
acts as an "endogenous" forcing able to permanently sustain the motion even in presence of 
friction. 

The previous interpretation of Lorenz dynamics actually leads to a generalization of the 
Lorenz model to a wider class of systems showing similar dynamical properties. We shall 
introduce a particular system belonging to such class which, because of its simplicity (piece- 
wise linearity), will allow us to study the dynamics of the original model using analytical 
tools. Indeed, in the steady chaotic regime (i.e. when the system permanently lies on its 
attractor set) the evolution of a point in the Lorenz phase-space can be roughly described as 
a somewhat regular amplified oscillation around one of its fixed points, followed by a sudden 
entry in the control basin of a twin fixed point symmetrical with respect to the origin, by a 
new amplified oscillation around the latter point, by a sudden return in the control domain 
of the former, and so on. 

The most evident aspect of the chaotic regime is just the unpredictability of the instant 
at which the center of the aforementioned amplified oscillations changes. The choice of a 
piecewise linearized version of the original model will allow us to highlight this point, while 
keeping unchanged the peculiar topological properties of the Lorenz dynamics. The exact 
equation ruling the instant of change of the oscillation center will be derived and a discon- 
tinuous dependence of this instant on the initial conditions will be highlighted. It will also 
be possible to write the analytical equations that defines the first-return two-dimensional 
Poincare map for the piecewise linearized system, which in turn synthesizes the main chaotic 
features of the model dynamics. Moreover, we shall show that, under suitable and reason- 
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able conditions, the evolution of the system completely reduces to a one-dimensional chaotic 
map. 

To summarize, starting from our interpretation, we shall be able to propose a piecewise lin- 
earized version of the Lorenz model, which on one hand has the same dynamical properties 
of the original system and, on the other hand, will provide analytical tools to explicate the 
emergence of chaos in Lorenz-like systems. 



II. THE LORENZ EQUATION 



The original Lorenz system [|l| consists of the three first order ordinary differential equa- 
tions 

X = -aX + aY 

Y = -Y + {r- Z)X (1) 
Z = -bZ + XY, 

where the dots indicate time derivatives and a, b and r are positive parameters originally 
related to the fluid properties and to the boundary conditions in the Rayleigh-Benard prob- 
lem. (Lorenz used a = 10, 6 = 8/3 and r = 28). 

The fixed points of system and their corresponding stabilities depend on r. For r < 1 
there is only one (stable) fixed point in (0, 0, 0). For r > 1 the origin looses its stability and 
a pair of new fixed points appear: (±[b(r — l)]"*^^^, ±[6(r — 1)]"*^^^, r — 1), which are stable in 
the range 1 < r < Vc with Vc = a{a + b + 3)/{a — b — 1). (Note that the critical value exists 
if and only if cr > b + 1). For r > all fixed points are unstable, and the Lorenz system can 



exhibit either periodic or chaotic behavior (see, e.g., Ref. ||2^ for a comprehensive exposition 
on the matter). 

Since we are interested in the case r > 1, we can suitably define the scaled coordinates 

X = X/[b{r — 1)]^/^, y = Y/[b{r — 1)]""^/^, z = Z/{r — 1), so that the system (|l]) becomes 

X = (y{y — x) 

y = —y + x + (r — 1)(1 — z)x (2) 
z = b {xy — z) 

with fixed points (0, 0, 0) and (±1, ±1, 1). 

We now reduce the Lorenz system to a unique differential equation for x = x(t), whose 
solution makes it possible the direct calculation of y{t) and z{t). By inserting into the third 
equation of (0) the expression for y in terms of x obtained from the first, one easily gets 



z + bz = — 
2a 



dt' 



— (x^) + 2ax^ 



(3) 



whose general solution is given by 



z{t) = e"''^*-*") ( z{to) - ^AU)) + ^^\t) + (l - ^ ) W d-^ e-'^'-^^x\s). (4) 



to 
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We have already stated that we are interested in the evolution of the Lorenz system in the 
chaotic steady state, i.e. far from the initial transient. Thus, we let to move back to — oo 
and obtain the steady-state expression for z{t), which can be written in the form 

where we use the square brackets notation [f]k{t) to indicate the exponential average of any 
suitable time function / on its past history: 

POO 

[Mt)^k dse-'^fit-s). (6) 
Jo 

From now on we shall refer to [f]k as the k- exponentially vanishing memory of the function 
/ or in short its memory. Note that, for the memory [f]k to exist, it suffices that f{t) ~ 
0{exp{—ht)), with h < k, a.s t — > — oo (the time functions we shall be dealing with are 
even bounded in this limit). 

To summarize, we insert in the second equation of (0) y{t) obtained from the first equation 
and z{t) given by and get for the variable x = x{t) alone in chaotic steady state the 
following second order differential equation with memory (in fact, an integral-differential 
equation) 



X + {a + l)x + cr(r — 




One can interpret the previous equation as the dynamical equation of a (unit mass) 
particle, viscously moving in a compound potential energy field consisting of a weighted 
average (in the chaotic regime b/2a < 1) of a quartic potential energy field independent of 
time (Fig. 0) 

U{x) = a{r-1)^-^^ (8) 
and of a quadratic potential energy field 

U,{x) = a{r-l)[x'-l],^^^, (9) 

whose curvature is given by a suitable memory function of the past motion. 

Without the memory term the particle would stop in one of the minima of the bistable 
potential U due to the damping term —{a+ l)x. Indeed, the Ut contribution yields, through 
an exponential average on the past evolution, an "endogenous" forcing term which can 
permanently sustain the motion. The particle oscillates with growing amplitude around the 
minimum of one of the potential well until its energy is sufficient to allow the crossing of 
the barrier in x = 0. As already said, the chaotic behavior of the system emerges just in 
the unpredictability of the instant at which the particle moves from one well to the other. 
Trajectories relative to very slightly different initial conditions can produce strongly different 
sequences in the number of oscillations in each well. In the next section we shall explain 
this fact analytically. 
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In order to simplify and to make the analysis standard, let us rescale the time t in Eq. 
as t — s> r = [(r — 1)6/2]^/^ t, to obtain the equation (hereafter called Lorenz equation) 

3^ oc doc 

-^ + Tlj^ + {x^ -l)x = -a[x^ -l]px, (10) 



where 



v/(r- 1)6/2' h ' ^ Vr-1' 

Note that, given h and r, a (the viscosity parameter in the original problem) affects both 
the friction term and the forcing term in the Lorenz equation. This fact shows how much 
these "opposite" contributions are in fact strictly related if the equation is to be viewed as 
representative of the original Lorenz system. Even if one now considers the Lorenz equation 
as the main subject of the study, one must note that not all the (positive) values of a, /3, r\ 
are consistent with their definition in terms of the original parameters 6, cr, r. In particular, 
in order to observe chaotic behavior, the following inequality must hold 

2/? ^ ' 

(Further details on the relation between the two sets of parameters are given in Appendix 
0). 

Equation (|10|) allows us to highlight the role of the memory forcing term in the Lorenz 
system dynamics. For this purpose it is interesting to compare the Lorenz equation with 
other examples of chaotic nonlinear (nonautonomous) systems as, for instance, the inverse 
Duffing equation x+?7x + (x^ — l)x = Acos(f2t), which describes a sinusoidally forced quartic 
oscillator or more appropriately the parametrically forced equation 



X + r/x + (x — l)x = — A cos(fit) X. (12) 

Note that in the latter cases the motion is sustained by externally assigned forcing terms, 
while in the case of the Lorenz equation the motion is self- sustained by the endogenous term 
—a [x^ — 1]^ X. In Figs. |^ and ^ we give a numerically obtained comparison between the 
phase portraits of the Lorenz equation and of equation (|1^) with parameters (a, /?, ?]), and 
(A, f2) suitably chosen in order to get similar ranges of motion. Note that the endogenous 
Lorenz forcing term, mimed by — Acos(fit) x, is in fact neither monochromatic nor with 
vanishing average value. 

III. THE GENERALIZED LORENZ SYSTEM 

The Lorenz equation (ITO) can be usefully recast in the form 



c? X (ix / \ 

^ + r/— + (g(x) + a [g(x)]^j <l>'(x) = (13) 

where the prime indicates the derivative with respect to x and in our case $(x) = 1/2 x^ 
and g(x) = x^ — 1. Such an equation can be interpreted as the description of the motion of 
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a unit mass particle subjected to a viscous force —rjdx/dT and interacting with a potential 
field $ through a "dynamically varying charge" Qrix) = q{x) + a[g(x)]/3. This charge both 
depends on the instantaneous particle position x{t) (by means of the term q{x)) and on 
the past evolution (by means of the memory charge [g(x)]/3). The coupling of [g]^ with the 
fixed potential field $ acts as an endogenous forcing term which can sustain the motion 
even in the presence of friction, and the chaotic behavior can actually arise from the synergy 
between this term and the viscosity. 

Put in the form (|13]), the Lorenz equation is arranged to be generalized to a generic charge 
q{x) interacting with a generic potential field $(x). Correspondingly, it is possible to obtain 
a generalized Lorenz system whose "x - projection" (far from the initial transient) yields 
Eq. ([13|) . Indeed, by inverting the calculation followed to derive Eq. ([T3| ) from system @, 
one easily gets the generalized Lorenz dynamical system 

' X = cr(y — x) 

< y = -y + X + {r - - z)^'{x) (14) 

^ z = -bz + b[^q'{x){y - x) + q{x) + 1]. 

Therefore, the specific Lorenz model can be viewed as singled out from a quite general class 
of dynamical systems which can exhibit chaotic behavior, their common essential property 
being an exponentially vanishing memory effect together with a viscous damping. 

Equations (|13D and ([T^) are related to Eq. (^ by assuming as potential energy field the 
quantity U = {b/2a)U + (1 - b/2a)Ut with 

U = a{r ~1) J q{x)^'{x)dx and Ut = a{r - l)[q]f3 ^ . (15) 

Obviously, any choice of q and $ should maintain the main properties of the Lorenz model, 
i.e., correspond to a two-well piecewise differentiable potential energy U{x), such that 
U{x) — > cxD as |a;| oo. 

We shall now focus our attention on a particular choice for q and $, that will maintain 
all the qualitative properties of the Lorenz system and, at the same time, will allow us to 
deal with chaos analytically. 



IV. THE PIECEWISE LINEARIZED LORENZ SYSTEM 

A. Linearization near fixed points 

As already noticed, the chaotic behavior of the Lorenz system essentially depends on 
the unpredictability of the instants when x change its sign: as long as it keeps constant 
sign the system evolution is certainly nonlinear, and nevertheless not "chaotic" at all. This 
fact suggests a slight modification of the original form of the Lorenz system, in order to 
single out analytically the origin of chaos without to be faced with the difficulties arising 
from nonlinear features. We thus set in Eqs. ([131) and (|T^) $(x) = |x| and q{x) = \x\ — 1 



obtaining (for x ^ 0) the piecewise linear Lorenz-like equation 



d'^x dx 

+ V-r + 'I ~ 1 + [ ~ 1 ]/3 I" sgn(x) = 0, (16) 
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where sgn(a;) = |a;|/x. The corresponding piecewise hnearized dynamical system, with the 
original choice of parameters, is then given by 

' X = a{y — x) 
y = -y + X + {r - - z) sgn{x) 

z = —oz + osgn(x) — - — . 

Our assumptions on $ and q correspond in Eq. (|^) to U{x) = a{r — 1) ( |a;| — 1 )^ /2 and 
Ut{x) = a{r — 1) [ |a;| — 1 ]^ ( |a;| — 1 ). Thus, we are faced with a simplified model, obtained 
by replacing the constant-in-time quartic potential with a piecewise quadratic one resulting 
from the superposition of two parabolas with vertex in ±1 and truncated at x = (Fig. |[). 
The two-well character of U is obviously maintained as well as the piecewise differentiability. 
The replacement of the original potential actually corresponds to a linearization of the 
system around both unstable fixed points (±1, ±1, 1), with the metching performed in a; = 
0. It appears that the chaotic behavior of the original model does not depend on the 
differentiability in x = 0. One can guess that also other classes of dynamical systems can 
be transformed in a piecewise linearized version by means of the same operations. 

It is easy to check that the fixed points of (0) are (±1,±1,1) and that the equation 
ruling their local stability is the same as for the original Lorenz system with parameters b, 
a, p = {r + l)/2. In particular, if a > 6 -|- 1, the critical value of r for our piecewise linear 
system ([T7| ) is given by rc'*"'* = 2rc — 1. In Figs. ^ and ^ two chaotic phase portrait for 
systems and (|l3) are shown, corresponding to the same choice of b and cr, and different 
choices of r in order to preserve the relationships between r and the proper corresponding 
critical value. 

When dealing with system (p!7|), the main simplification is that it is separately linear for 
X < and x > and can thus be analytically solved in each region. Indeed, by applying 
the operator {d/dr + (3) to each side of Eq. (|TB[) one obtains the equation 

^ + (/? + + (1 + (3v)^ + /3{l + a)ix- sgn(x)) = , (18) 

which can be explicitly solved separately on each side of x = 0. The nonlinearity of the 
model is simply reduced to a change of sign of the forcing term ±/3(l + a) when x crosses 



the plane x = 0, henceforth denoted with vr As we shall see, the crossing times are 



somewhat unpredictable, as they result from the discontinuous solutions of an (incidentally 
transcendental) equation. Our piecewise linearized system will thus turn out to be an 
important tool to analytically investigate the emergence of chaos in Lorenz-like systems. 
The main advantage of the piecewise linearization is that one has to deal only with the 
simplest nonlinearity, i.e. an isolated (not eliminable) discontinuity (see, e.g., Refs. [P^-p7| 
for other examples of piecewise linear chaotic dynamical systems). 



B. Analysis of the motion 

Let us now consider in detail the second order integral-differential equation (p!6D, which 
describes the evolution in time of x for the piecewise linearized system in the chaotic steady 
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state. It is equivalent to a third order nonlinear differential equation whose phase space is 
described by the coordinates x, x, x (with a small abuse of notation from this time on we 
shall indicate with the dot the differentiation with respect to r). Solutions of Eq. ( [T6| ) can 
be easily calculated with tools of the customary analysis on each side of x = 0. To obtain 
a global solution, such partial solutions should be matched at a; = under the reasonable 



assumptions that the position x, the velocity x, and the memory [ |x| ]/3 are continuous p8 
Notice that, in contrast, when crossing the plane tt, the acceleration x turns out to be 
undefined. However, if x{f) = 0, referring to the left and right time limits x{f~) and x{f~^), 
it appears from Eq. (|T^ that they are related by the equation x{f~) + x{f~^) = —2rjx{f). 
In the sequel we shall refer to x{f~) and x(f+) as the acceleration "immediately" before 
and "immediately" after the crossing, respectively. 

Since Eq. (^) is invariant under the transformation x — > — z, r — > r, we can focus our 
attention only on one of the two regions, e.g. a: > 0, and describe the motion in this half- 
space (the evolution in its twin half-space being recovered through the change x —x). 
Clearly, this is equivalent to put a "rigid wall" in x = and look at the crossing of vr as an 
elastic collision. 

In summary, according to the previous scheme, the system evolution in time is completely 
described by the following steps: 

1. motion for a; > 0; 

2. collision against vr and discontinuity of x; 

3. inversion x —x and matching with a new solution defined in the region x > again. 



1. Motion in the half-space x > 
If we define ^ = x — 1 , Eq. (|T8|) assume a simple form 

+ + + + /9(l + «)e = 0, (19) 

that is a linear third order differential equation (homogeneous and with constant coefficients). 
As a consequences of the Routh-Hurwitz theorem , a critical value ac for the parameter 
a exists, i.e. 

a„=(i±^|^-l. (20) 

For a > ac the fixed points are unstable saddle focus with a real negative eigenvalue — Aq 
and a pair of complex conjugates eigenvalues Ai = A,. + Ai = A,. — iXi. 

For the sake of simplicity we shall indicate with (0, ri) the time interval between two 
consecutive collisions in the steady-state of the system. Without loss of generality we assume 
X = and x > 0. Then, it can be easily shown that in the time interval (0, ri) the motion 
in the phase space ^, C,, > —1) is completely described by the equation 

^(r) = M{t) M{0)-' , (21) 
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where we have defined 






/ -1 \ 




m 







and 



-Aqt 



M(r) 



( Re(e^i^) Im(e^i^) 

Re(Aie^i^) Im (Ai e^i^) -Aq e"^"^ 
\Re {XI e^i") Im (A? e^^") Ag e"^"" ) 



(Note that, for a > M.{Q) is always invertible since det A^(0) = Aj [Af + (A^ + Aq)^] and 
Aj > 0). The eigenvalues of the matrix M.{t) A^(0)~^, which connects the vector ^(r) to its 
initial value , are e'^^^, e^^'^ , e~^°'^, with corresponding (constant-in-time) eigenvectors 



Vi 













Ai 










With respect to the base {vi, V2, v^} one has 

^(r) = e^^^ (ci e^^^vi + C2 e-^'^V2) + C3 e'^'^v^ . 



(22) 



For the saddle point ^ = there exist a stable one-dimensional manifold corresponding 

to V3 

W = {pvs I p > -1} 
and an unstable two-dimensional manifold generated by Vi, V2 (see Fig. ^ 



Vi + V2 



Vl - V2 

2i 



p > —1 , g G 



The evolution of the linearized system in the interval (0, ti) is the combination of an ex- 
ponential decay along and of an amplified rotation on W". To determine the relative 
quickness of each component of the motion with respect to the other, consider that from 
Eq. (|l^) it is easily checked that Aq — 2Ar. = P + 1], and so Aq > A^ V a, P, rj. Therefore, the 
exponential decay along the stable manifold is always more rapid than the exponential 
oscillating growth on the unstable manifold W", the parameter controlling this difference 
being jS + rj. 

Although the fixed points are saddle focus, it appears that the piecewise linearized system 
does not exhibit Shilnikov chaos owing to the absence of an homoclinic orbit. 
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2. TT-coHision and inversion 



To complete the description of the piecewise hnearized system dynamics, we consider 
the instant ti at which the first colhsion occurs. At this time, as aheady remarked, the 
trajectory coming from the half-space a; > must be matched with the solution defined 



again in the same region, but with corresponding new "initial conditions" [31 



or, in matrix form. 



with 



-1 

-iin) (23a) 
-e(ri+) = |(rr) + 2r^e(ri), 



^1 = ai?^(r+) (23b) 




V=\01 01 and 

The matrix V accounts for the acceleration discontinuity in x = 0, while 3 yields the 
sign inversion after the impact. From Eqs. (|2l|) and (p3b|) it follows that the velocity and 
the acceleration immediately after the collision are related to the initial conditions by the 
operator V{t) = JV Ai^r) A4{0)~^, according to the formula 

^i = ^(n(O)^0- (24) 

Notice the highlighted dependence of ti on ^q, which reveals the nonlinear character of this 
important relationship. 

Starting from the above results, the very origin of chaos in the piecewise linearized Lorenz 
system will be identified and discussed in the next section. We shall also show that the basic 
mechanisms for chaos to emerge apply also to the original Lorenz system. 



C. Dependence of ri on initial conditions 

1. Unpredictability of the crossing time 

As already observed, the instant at which the plane n is crossed, strongly depends on the 
initial conditions. This fact is strictly related to the chaotic behavior of Lorenz-like systems. 
Let us now study in some detail this topic for the piecewise linearized model. 

The instant ri at which the first collision occurs is defined by the condition ^{ti) = — 1. 
From the first line of Eq. ( pl| ) ti is thus the smallest positive solution of the transcendental 
equation 

-1 = Ci e^^^' cos(Airi) + C2 e^^^' sin(Airi) + C3 e"^""^ (25) 
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where Ci, C2, C3 are linearly related to initial conditions. The residence time ti is therefore 
the first intersection of the graphs of giji) = —C3 e~^°'^^ — 1 and h{Ti) = Ci e'*"''^^ cos(Ajri) + 
C2 e^^'^^ sin(Ajri). Since is a decreasing exponential function and h an oscillating function 
with growing amplitude, one can easily understand why even a little modification of initial 
conditions can produce a discontinuous variation of ri (Fig. ||). 

In conclusion, the unpredictability of the residence time is closely connected to the dis- 
continuous character of the solutions of Eq. (pS]). The chaotic behavior for the the piecewise 
linearized system clearly stems from such unpredictability. Our claim is that an entirely 
analogous situation exists also for the original Lorenz model. We shall return however in 
more detail on this important analogy in the sequel. 



2. The residence time ti as a function of ^0 <ind ^0 

We now explicitly investigate the dependence of ti on the initial velocity and acceleration. 
As previously remarked, the function ti = ti{C,o, ^0) is implicitly defined as the smallest 
positive solution of Eq. (^^). Since Ci, C2, C3 are linearly related to the initial conditions, 
Eq. (p5D can be rewritten in terms of ^0 and ^o- In this form it describes a family S of 
straight lines in the plane {^o;^o} parametrized by ti 

S : A{n)io + B{n)^o + Cin) = (26) 

[see Appendix for the explicit expression of the coefficients A{ti), -B(ti), C{ti)]. 
Let us denote with (Ti), i >0 the ordered sequence of the zeros of B{ti). Note that one 
always has Tq = 0. The slope of the straight lines of 5, —A/B, and their ordinate for ^0 = 0, 
i.e. —C/B, have their singularities in Tj, i > 0. Both these functions are only asintotically 
periodic because of the presence of terms proportional to e~'^'''^\ which become negligible 
only for large Ti (Figs. ^ and pUj ). 



In each interval (Tj, Tj+i), i > 0, the function —A/B is everywhere growing (see Appendix |B[) 
and varies from —00 and +00. Therefore, as ti increases, the straight lines of S rotate 
anticlockwise and at the same time translate starting from the ^o-axis (obtained for ri = 0) 
(Fig.0). 

The envelope of the family iS is a curve 




which looks similar to an elliptic spiral (Fig. 12). (A parametrical representation of 7 is 
given in Appendix P) . 

From Eq. (^) the contour lines of the function ti are parts of straight lines in the plane 
(^0; ^o)- Indeed, according to the definition, each line of S is in fact a contour line for ti only 
where its points do not belong to another line corresponding to a smaller value of ti and, 
furthermore, only where their abscissas correspond to positive values of .^o- Rather than a 
family of straight lines, the function contour lines are raies or segments according to the 
constant value of ri. 

For Ti ranging from to T2 these lines perform a complete anticlockwise "rotation," starting 
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from the ^o-axis. Thus, they cover the entire half-plane > except the region inside their 
envelope; from this first rotation one obtains a set of raies (Fig. O). 

Subsequently, for ti ranging from T2 to T4 the contour lines become segments of variable ori- 
entation and contained in the region delimited by the curve already generated from the first 
rotation and the new envelope of the second set of straight lines. This behavior repeats itself 
V Ti G [T2i,T2i+i], i > 0. Because of the aperiodic character of the functions involved, the 
curve 7 outlines in the positive half-plane a structure consisting of "pseudo-elliptic" annula. 
Moving anticlockwise along each of these annula, ti grows continuously. On the contrary, 
passing through the border that separates two different bands, one meets discontinuities in 
the dependence of ri on initial conditions. 

The curve 7 "winds" round the point Pq = (Aq, — Aq) (see Fig. |T^. For {C,o,^o) (Ao, — Aq) 
one has ri ^ 00, since for these initial conditions one obtains Ci = C2 = and the system 
exactly lies on the stable manifold W*. Its motion is in this case an exponential decay 
towards the fixed point ^ = 0. 

The previous observations allow to easily guess the structure of the graph of ri = ri(^o, ^0) 
shown in Fig. It should be noted that ti shows instability with respect to the initial 
conditions only in a limited subset of the half-plane ^0 > 0. As we shall see, it is natural to 
expect that in the chaotic regime the system is quickly attracted inside this region. 



D. The piecewise linearized model as a one-dimensional map 

1. The TT -plane Poincare map 

At this point of our study we have analyzed in some detail the unpredictability of the 
time at which the system crosses the plane vr. To completely motivate the chaotic dynamics 
of the piecewise linearized system and thus of the original model, we must, however, add 
some further results on the attracting set of the system. 

All results obtained in the previous section can be easily extended to the n-th colli- 
sion against vr. Specifically, denoting with ^„ the array assigning position, velocity, and 
acceleration immediately after the n-th collision 




we have [similarly to Eq. (|2^)] 



€n=nrn)^„-l, (27) 



where we have denoted with r„ the n-th residence time, i.e., the time interval between the 
{n — l)-th collision and the n-th one. 

The first line of Eq. ( pTf ) gives r„ in terms of ^n-i and C,n-i- the dependence of this time 
on the {n — l)-th initial conditions has already been discussed in Section IVC2| . Knowing 



r„, the other two lines allow to relate the velocity and the acceleration ^„ to and 
C,n-i ■ Eq. ( p7D defines a two-dimensional first-return Poincare map between (^„,^„) and 
(^„_i, .^„_i), obtained from the section ^ = — 1 of the phase space ^, ^, C,. We recall that, in 
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spite of its appearance, the map (|27|) is in fact non-linear, since the matrix V depends on 
r„, which is a transcendental discontinuous function of the (n — l)-th initial conditions. 
Indeed, under suitable conditions, the two-dimensional map for the linearized system prac- 
tically reduces to a one- dimensional map. 



2. Approximation of the system with a one- dimensional map 
It has been already shown that for ^ > — 1 the motion in the phase space has the 



form dl) 



where Aq and A^ = Re(Ai) satisfy the equation Aq — 2Ar = P + t], which in turn implies 
Ao > Xr- As a consequence the evolution in the phase space ^, ^, ^ consists both of a "rapid" 
decay towards ^ = along the stable manifold and of a "slow" amplified oscillation on 
W". We thus expect the phase trajectories to be strongly attracted on the unstable manifold 
and, once on W", to slowly spiral outwards. Provided trajectories start close enough to W", 
then they meet the vr-plane very close to its intersection with W" itself, i.e., along the 
straight line 

+ 2« + (A; + A?) = 

Let us assume that the attraction towards the unstable manifold is very strong and thus it 
takes place almost instantaneously (the goodness of this assumption is controlled by l3 -\-rj). 
Under this hypothesis all trajectories approximately hit vr along the straight line C~. Thus, 
from Eq. ( p3b|) , it follows that, immediately after each collision, the system necessary lies 
very close to the straight line 

|e + 2(r7 + A.)^-(A2 + A?)=0 ^^g) 

Therefore, the following relation between velocity and acceleration is expected to approxi- 
mately hold 

4 + 2{7] + Xr) in - (Xl + A^) = . (30) 

From Eqs. (^) and ( PDD an analogous linear dependence between the velocity and the 
memory follows 

iv + 2Xr)L -aWn + l-{Xl + X^) = , (31) 

where Wn denotes the /3-memory of ^ evaluated at the n-th collision. As a consequence, in 
the limit we have considered, the attracting set for the map (p7|) is C'^ and, therefore, it 
reduces to an one-dimensional map. 

Since in the chaotic steady-state the velocity and the acceleration that define the initial 
conditions for the trajectory after each impact are not independent, the time can be 
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expressed as a function of ^„_i alone: = r„(,^„_i) (see Fig. |T^). The behavior of this map 
is easily understood if we refer to the graph of ri as a functions of and i^o (see Fig. 
Indeed, in the chaotic regime is superimposed to the region of the plane (^„_i,^„_i) 
contained by 7, where ti shows unstable behavior with respect to a change in the initial 
conditions (Fig. |16D. For those values of ^n-i, for which passes through the same "pseudo- 
elliptic" corona, r„ slowly changes with varying crossing velocity. On the contrary, when 
intersects the boundary between two different coronas, r„ shows a discontinuity in its 
dependence on (Fig. plsl) . 

The linear dependence between velocity and acceleration immediately after each collision 
implies that the Poincare map ( ^Tf ) becomes one-dimensional, e.g. a map between ^„ and 
Figure ^ shows this fact for a = 6.50, (3 = 0.19, r] = 0.78. By simple inspection of 
the derivatives corresponding to the fixed points of the map one can easily check that they 
are all unstable, so that the map produces a chaotic behavior. Note that the discontinuities 
simply correspond to the analogous ones in the function = r„(^„_i), evaluated along the 



straight line given by (pOD and drawn in Fig.|16 



For the sake of completeness, we turn now to the description of the piecewise linearized 
system in the original phase space x, x, x. There the system has two fixed points (±1, 0, 0). 
For each one there exists a stable manifold 

W'^ = {pvs+ (±1,0,0) I p^Tl} (32) 

and an unstable manifold 

>V^ = + ^^Sr^ + ' p^TlandgGMj . (33) 

During the evolution the trajectories are rapidly attracted on the unstable manifold relative 
to the half-plane where they belong (Fig. |T8|). Thus, immediately after the crossing of the 
TT-plane the phase point lies on one of the straight lines 

'x + 2{r] + Xr)x + (A^ + A^) sgn(x) = 
X = 0, 

while the memory and the velocity are linearly related according to equation 

{r] + 2A,) Xn - a sgn(i;„) «;„ + 1 - (A^ + A^) = . (34) 

As for the crossing of the vr-plane, the system can therefore be described by a one-dimensional 
map. To get the maps Tn = rn{xn-i) and i;„ = x„(i;„_i) from the analogous ones in the 
rigid wall scheme, one has simply to consider that ^„ coincides with \xn\- 

The reader should note that all the maps we have drawn in the previous figures have 
been obtained by numerically solving the analytical equation (pTl), which exactly defines r„, 
and ^n+i, and they have not been computed by numerical integration of the differential 
system (as usual for non solvable dynamical systems like the Lorenz original one). Thus, 
the use of numerical tools has been required only because of the transcendental character of 
the concerned equation. 
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We explicitly remark that in general the attraction towards is not infinite, but 
nevertheless more rapid than the amplified rotation on the unstable manifold (remember 
that Ao — 2Ar = P + ri>0\/a,j3,ri). Therefore, we can expect that, if /3 + 77 is finite, the 
attractor set is not exactly the straight line C~^, but a narrow "strip" that contains jC~^. This 
strip intersects the region of the graph of r„ where the residence time is strongly dependent on 
initial conditions and the mechanism for chaos to arise is absolutely the same. Our numerical 
simulations for finite (3 + 7] are in excellent agreement with these predictions and then we do 
not show their corresponding graphs here, since they are practically indistinguishable from 
the graphs we have already discussed. 



V. COMPARISON WITH THE ORIGINAL LORENZ SYSTEM 

We conclude our analysis by showing some simulated numerical results. Bfore doing 
that, some remarks are worth discussing. We previously considered independent 
coordinates in the phase space because of their physical meaning: it is more intuitive to speak 
about acceleration rather than of memory of the system. Unfortunately, if we want to "assign 
initial conditions" to the original Lorenz system after the crossing of the plane vr, we cannot 
consider the acceleration and the velocity since from Eq. ([10|) they are not independent 
variables in x = 0. Indeed the relation Xn + TjXn = always holds, while the memory is 
not known a priori. As a consequence, we have to consider as independent initial conditions 
position, velocity, and memory. However, it is worth remarking that the comparison between 
the Lorenz system and the piecewise linearized one is completely meaningful. This because 
for the latter model Wn-i is easily known in terms of and from Eq. (T^): for the 
linearized system the choice of C,n-i or of Wn-i is absolutely equivalent. 

Since we are interested in the behavior of the Lorenz system in the chaotic steady state, 
all the following results are thus obtained via the numerical integration of the steady-state 
differential system 

' X = V 

V + rjv + (x^ — l)x = —a w x (35) 
w = -Pw + Pix"^ - 1), 

where now w = [x'^ — 1]/?. This dynamical system, as already proved, is equivalent to the 
system (|^) far from the initial transient. 

We can now turn to the comparison with the piecewise linearized system. For the original 
Lorenz model we have computed the maps that relate and \xn\ to |xn-i|, respectively, 
by numerical integration of the system ( pS]) (the notations are the same as in the previous 



section). Figures |T9| and |2y show the results we obtained: the qualitative behavior of these 
maps is very similar to that of the analogous maps for the piecewise linearized system. As 
we will see, the mechanism for chaos to arise is indeed the same. 



Observe that the maps in Figs. |T9] and |2y are not perfectly onto: obviously this is due to 
the fact that, in general, the Lorenz system does not reduce exactly to an one-dimensional 
map when the plane a; = is crossed. 

From direct computation it appears that the residence time around a fixed point depends 
on the initial conditions in a way similar to that of the linearized system (Fig. |2l|). There is 
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a region in which ri varies continuously with varying initial conditions and a region (Figs. ^ 
and |2^) of strong instability. In the latter region ti changes continuously along each "annu- 
lus" of the "spiral," whereas it shows abrupt discontinuities crossing the boundary of each 
annulus. As expected, the system is attracted on that region (Fig. |2^) and the maps in 
Figs. O and EOl derive from the overplotting between Figs. E3] and |22|. 



The chaos actually emerges due to the combination of the step-like first-exit-time (Fig. jTOj) 
and of the piecewise return map (Fig. ^U[) for initial conditions. This mechanism is the 
same as for the linearized model, for which we gave an analytical description. The results 
we obtained for the piecewise linearized system therefore provide in our opinion a useful 
enlightening tool to understand both the original Lorenz system and its eventual general- 
izations. 

Finally, in order to point out the connection between the original Lorenz system and its 
piecewise linearized version, it is illuminating to compare the tent-like map shown by Lorenz 
in his original work with the analogous map for our simplified system. Describing the 
relation between two consecutive maxima of the coordinate z{t), Lorenz found a univocal 
tent-like map (Fig. [2^ ). For our linearized system we obtain similar results (see Fig. p5[ ). 
The growing branch of the map is approximately linear, since for small values of z the system 
evolution is confined in a half-space and, there, it essentially consists of an exponentially 
amplified oscillation. For large values of z the map has the same shape of the Lorenz tent- 
like map. The linear growing of the left branch of the map is due to the simplification we 
made on the evolution on each side of x = 0. The shape of the right branch shows that 
the piecewise linear system keeps all the complexity of the Lorenz model, which can be 
attributed to the unpredictability of the x = crossing, which is the very origin of chaos in 
Lorenz-like systems. 



VI. CONCLUSIONS 

The main conclusions from the present paper can be summarized as follows. We have 
reformulated in a classical mechanics form the dynamics of Lorenz-like systems and showed 
that its three-dimensional phase-space dynamics can be mapped into an one-dimensional 
motion of a particle oscillating in a conservative quartic two-well potential, subjected to 
a viscous dissipation and to a memory forcing. Starting from this interpretation, we have 
introduced a piecewise linearized version of the Lorenz system (belonging to a larger family 
of Lorenz-like systems) which substantially has the same properties of the original model 
with the advantage that it allows an analytical treatment. 

The most evident aspect of chaotic regime is the unpredictability of the instant at which 
the center of the particle amplified oscillation changes. Chaos arises due to the combination 
of the step-like behavior of this time with the piecewise return map defining the crossing 
conditions. This aspect has been singled out analytically, focusing on the piecewise linearized 
version. There, the exact equation for the time, at which the oscillation center changes, has 
been derived and the discontinuous dependence of this time on the crossing conditions has 
been shown analytically. 

By means of numerical simulations we have verified that the highlighted mechanisms for 
the chaos to emerge survive also for the fully-nonlinear Lorenz system, where analytical 
techniques are not applicable. 
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APPENDIX A: PARAMETERS a, /3, r/ 



In the Lorenz system, the following conditions for cr, b, r have to be satisfied in order to 
have three instable fixed points: 



a, 6, r G M 
r > 1, 
a >b + l, 
r > Tc . 

The parameters a, (3, rj are defined as follows: 

2a 



(Ala) 
(Alb) 
(Ale) 
(Aid) 



a = --l, 



2b 



r — 1 



V 



CT+1 



^/{r-l)b/2' 



From (|Ala|) and ( |Alb|) it immediately follows that a, f3, rj are positive too. Therefore, a, 6, r 
can be rewritten in terms of a, /3, 77 in the form 

(l + a)/3 



r = 1 



/3(a + l) -2r/' 
4 

/5(/3(a + l)-2r/)' 
/3(a + l) -2r/' 

Under the hypothesis a, P, r] E R^, ( |Ala| ) and (|Alb|) are equivalent to 

Pia + 1) < 2ri, 

while (|Alc|) corresponds to the inequality 

(3{a + l) > (3 + 7]. 

Specifically, (|Ala| ), ( |Alb|) , and (|Alc|) then imply 

V 

a>->l. 



(A2) 



(A3) 



(A4) 
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Finally, ( |Alcl|) can be rewritten in terms of new parameters as 



a > 



2/3 



(A5) 



from which necessarily ([A4|) follows. 



APPENDIX B: FAMILY S 

The family of straight lines S is defined by the equation 

S: A{n)io + B{n)^o + C{n) = 

where 

A{ti) = -2AiA,e-^o^i + e^'-^i(2AiA, cos(Airi) + (Ag + A^ - A^) sin(Airi) ) 
B{n) = Xi e-^o-i + e^^^' ( (A^ + Aq) sin(Airi) - A^ cos(Airi) ) 
C{ti) = -A, (A^ + A^) e-^«-^ - Aq e^^^^ [ A, (2 A, + Aq) cos(A,ri) + 
+ (Af - (Ao + A,)A,) sin(A,ri) ] + ( A, (A^ + (Aq + A,)^ ). 

1. Slope of straight lines of 5 

We indicate with (Ti)i>o the ordered sequence of the zeros of -B : Ti < Tj+i 
V i > with B{Ti) = V i > 0. The first derivative of —A/B is positive over the whole 
domain of definition. Indeed one has 

_ d /A 
~d^i [b 



A.e^-"! [A^e^-'i + e-^°"^ ( - A, cos(A,ri) - (Aq + A,) sin(A,ri))] (A^ + (Aq + A 
A-e-^ori + gA.n [ _ cos(Airi) - (A^ + Aq) sin(Airi)]^ 

with 

oo 

Xie^^^' + e-^«"^ ( - Xi cos(Airi) - (Aq + A,) sin(Airi)) > V ri G [j{Ti, Ti+i) C 



i=0 



The latter inequality can be easily proved. Indeed: 

(A, + Ao)ri > log(l + (A, + Ao)ri) 
e(V+Ao)n > 1 + (A, + Ao)ri > cos(A,ri) + ^A±M 



A,: 



Xie^^^' > e-^°^' (A, cos(A,ri) + (A, + Aq) sin(Airi)) V n G |J(Ti, T,+i) C 



i=0 



2. Envelope of family <S 

A parametrical (not regular) representation of the envelope of the family S is obtained 
deriving and with respect to Ti. This can be done from the system 

r A{n)^o + B{n)^o + C{t,) = 



The result 
is given by 



and 



6(Tl) = 



7 : 



io = 6(ri) 
lo = 6(ri) 



{ - Ao(Ai(l + e-^°"0 cos(Airi) - Ai(e-(^°+^'-)^i + e^""!)) 
- (A^ + A,^ - AoA,)(e-^°-i - 1) sin(A,Ti)} ^ 

(AiC^""! + e-^°^i(-(Ao + A,) sin(Airi) - A^ cos(AiTi)) 



{A^(Ai(cos(AiTi) - e^'-^i) + A^sin(AiTi)) 
+ (e-^on _ i)(A.cos(A,ri) - A,sin(A,ri))(A^ + A,^) + 

Ao(e-^°^H(A2 - A^) sin(A,T) + 2A,A,cos(A,ti)) - 2A,A,e-(^°+^^)^i)} ^ 

(A^e^^^^ + e-^°^i(-(Ao + A,) sin(Airi) - A^ cos(A,ri)) . 
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FIGURES 



FIG. 1. The constant-in-time quartic potential f7 (cr = 10, 6 = 8/3, r = 28). The classical 
particle representing the Lorenz system moves in the potential U subjected to a viscous damping 
and to a memory forcing. The minima of the potential wells correspond to the unstable points of 
the three-dimensional Lorenz system. 

FIG. 2. Typical chaotic phase portrait for the original Lorenz model {r] = 1.31, a = 10.30, 
P = 0.216). 

FIG. 3. Typical chaotic phase portrait for the inverse parametrically forced Duffing equation 
(?? = 1.31, A = 4.2, = 0.99). The chaotic evolution of this system is characterized by the 
competition between a viscous friction and a forcing term with the same structure of the right 
hand side of the Lorenz equation. Here, however, the forcing is externally given and the memory 
is replaced by a known function of time. 

FIG. 4. Quartic potential U for the original Lorenz system (dashed line) and for the piecewise 
linearized system (full line) . The linearization of the Lorenz system maintains the qualitative shape 
of the constant-in-time potential. 



FIG. 5. Typical chaotic phase portrait for the original Lorenz system (a = 10, b = 8/3, r = 28). 

FIG. 6. Typical chaotic phase portrait for the piecewise linearized Lorenz system {a = 10, 
6 = 8/3, r = 55). 



FIG. 7. Stable and unstable manifolds for the point ^ = in the phases space ^, ^, ^. The 
evolution of the piecewise linearized system on each side of vr consists of an exponential decay along 
yV^ and of an amplified oscillations on W". 

FIG. 8. Graphical interpretation of the discontinuous character of ti for small changes of the 
initial condition. Full line is the curve g{Ti), dashed lines is /i(ti) for the parameters a = 6.50, 
P = 0.19, 7] = 0.78. Bullets denotes the first intersection between g and h, whose abscissa defines 
the residence time. 



FIG. 9. Slope of the straight lines of the family S {a = 6.50, /3 = 0.19, r/ = 0.78). [See Eq. (|2|) 
and Appendix p| for the exact definition.] 



FIG. 10. Ordinate for ^0 = of the straight lines of the family S {a = 6.50, (3 = 0.19, r/ = 0.78). 
[See Eq. ( [2^ ) and Appendix ^ for the exact definition.] 



FIG. 11. The family of curves S [defined in Eq. (|6|)] for a = 6.50, p = 0.19, and r] = 0.78. 
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FIG. 12. The envelope 7 of the family of straight lines S looks similar to a spiral {a = 6.50, 
P = 0.19, and r] = 0.78). 

FIG. 13. Contour lines of the residence time ri as a function of the initial conditions (^O; ■^o) 
for Ti G [0, T2] (a = 6.50, /3 = 0.19, r] = 0.78). 

FIG. 14. Graph of ti as a function of the initial conditions, obtained by numerical solution of 
the transcendental equation (p^) (a = 6.50, /3 = 0.19, r] = 0.78). 

FIG. 15. Map of the first-exit time function of the crossing velocity obtained by 

numerical solutions of Eqs. and (||). (a = 6.50, /? = 0.19, rj = 0.78) 

FIG. 16. During chaotic evolution the couple of crossing conditions {£,n,£,n) is attracted on the 
straight line , which is superimposed to the instability region of the residence time Tn (identified 
by the curve 7). The graph refers to the values a = 6.50, f3 = 0.19, r] = 0.78. 

FIG. 17. The velocity at the n-th crossing of the plane x = as a function of the velocity 
(a = 6.50, (3 = 0.19, t] = 0.78). Notice that, however, in the steady-state chaotic evolution 
the map is restricted to the range (0, 1.2) of possible values of S,n-i- 

FIG. 18. Attraction of the piecewise linearized system on the manifolds Wj.. The evolution of 
the system on each side of vr consists of a "rapid" exponential decay along the stable manifold and 
of a "slow " amplified oscillation on the unstable manifold. 

FIG. 19. Map Tn = Tn{\xn-i\) for the original Lorenz system, obtained by numerical integration 
of the system (||) (a = 10.30, p = 0.216, r] = 1.31). 

FIG. 20. Map between and [in-il for the original Lorenz system, obtained by numerical 
integration of the system (|35|) (a = 10.30, /3 = 0.216, rj = 1.31). 

FIG. 21. Graph of ri as a function of the crossing conditions for the original Lorenz system, 
obtained by numerical integration of the system ( |35| ) for different initial conditions (a = 10.30, 
P = 0.216, r] = 1.31), 

FIG. 22. Graph of ri as a function of the crossing conditions viewed from the above in the case 
of the original Lorenz system (a = 10.30, f3 = 0.216, r] = 1.31). 

FIG. 23. Map between the crossing memory Wn and the corresponding crossing velocity 
for the original Lorenz system (a = 10.30, /? = 0.216, r] = 1.31). 
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FIG. 24. The tent-like map for the original Lorenz system (b = 8/3, a = 10, r = 28). M, 
denotes the n-th maximum of the coordinate z{t). 

FIG. 25. The tent-like map for the piecewise linearized system (6 = 8/3, cr = 10, r = 100). 
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